Zitterbewegung of electrons in quantum wells and dots in presence of an in-plane 

magnetic field 
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We study the effect of an in-plane magnetic field on the zitterbewegung (ZB) of electrons in a 
semiconductor quantum well (QW) and in a quantum dot (QD) with the Rashba and Dresselhaus 
spin-orbit interactions. We obtain a general expression of the time-evolution of the position vector 
and current of the electron in a semiconductor quantum well. The amplitude of the oscillatory 
motion is directly related to the Berry connection in momentum space. We find that in presence of 
the magnetic field the ZB in a quantum well does not vanish when the strengths of the Rashba and 
Dresselhaus spin-orbit interactions are equal. The in-plane magnetic field helps to sustain the ZB in 
quantum wells even at low value of kod (where d is the width of the Gaussian wavepacket and fco is 
the initial wave vector). The trembling motion of an electron in a semiconductor quantum well with 
high Lande g-factor (e.g. InSb) sustains over a long time, even at low value of kod. Further, we study 
the ZB of an electron in quantum dots within the two sub-band model numerically. The trembling 
motion persists in time even when the magnetic field is absent as well as when the strengths of 
the SOI are equal. The ZB in quantum dots is due to the superposition of oscillatory motions 
corresponding to all possible differences of the energy eigenvalues of the system. This is an another 
example of mult i- frequency ZB phenomenon. 

PACS numbers: 75.70Tj,03.65.-w,73.21.La,73.21.Fg 



I. INTRODUCTION 

In recent years there is a growing interest in the field 
of spin based electronic devices. There has been a lot 
of study in this field after the proposal of the spin field 
effect transistor by Datta and Das-. The charge carriers 
carry spin in addition to their charges. The ultimate goal 
of this field is to control the spin degree of freedom of the 
charge carriers to produce and detect spin-polarized cur- 
rent in semiconductor nanostructures 2 -. One can develop 
device technology^ and quantum information processing 4 
in future on the basis of the manipulation of the spin 
degree of freedom. The coupling between intrinsic spin 
of an electron with its orbital angular momentum con- 
stitutes the intrinsic spin-orbit interaction (SOI) in low- 
dimensional semiconducting systems. Particularly, there 
are two kind of SOI present in low-dimensional semicon- 
ductor structures. One is the Rashba spin-orbit interac- 
tion (RSOI) and another is the Dresselhaus^ spin-orbit 
interaction (DSOI). The RSOI arises mainly from the in- 
version asymmetry of the confining potential in semicon- 
ductor hetero junctions. The strength of RSOI is propor- 
tional to the electric field externally applied which can 
be tuned by an external bias 7 -^ or internally generated 
due to the crystal potentials. On the other hand DSOI is 
present in bulk materials and semiconductor heterostruc- 
tures which lack bulk inversion symmetry. The form of 
DSOI term strongly depends on the growth direction of 
semiconductor quantum well (QW)&~— . The strength 
of DSOI depends on properties of the material and the 
crystal structure. 

In 1930, Erwin Schrodinger— predicted that a free 
particle described by the relativistic Dirac equation will 
perform an oscillatory motion which is known as the 



zitterbewegung (ZB). The free relativistic Dirac particle 
has two energy branches: e(p) = ±a/ to o 2 c 4 + c 2 p 2 . It is 
well understood that this oscillatory motion results from 
the interference between these two energy branches 1 ^. 
The large oscillation frequency lu c ~ 10 21 Hz and the 
small oscillation amplitude A c ~ 10~ 13 m is not accessi- 
ble to the modern experimental techniques. Most of the 
studies on the ZB of electrons used plane waves to de- 
scribe the electrons. It was pointed out by Lock-i 4 - that 
plane wave is not a localized state and therefore, rapid 
oscillations on the average position of a plane wave has 
some limitations. He also demonstrated that when an 
electron is described by a wave packet, the ZB oscillation 
has a transient character. 

In 2005, Zawadzki et alJ£ studied ZB m narrow gap 
semiconductor (NGS) by using the analogy between kp 
theory of the energy bands in NGS and the free Dirac rel- 
ativistic equation for electrons. They found much more 
favorable amplitude and frequency of the oscillation than 
those in a vacuum for a free electron. At the same time, 
Schliemann et aL 16 ' 17 has studied the ZB of an electron in 
III-V zincblende semiconductor QWs in the presence of 
the SOI thoroughly. The above studies initiated intense 
theoretical research on the ZB of electrons in various con- 
densed matter systems^— . 

The proposal of an experimental scheme for ob- 
serving ZB in ultra-cold^ atomic gases is also given. 
The trembling motion was proposed for photons in 
a two-dimensional photonic crystal and for Ramsey 
interferometryS 7 -. It was reported in Refi^S that an acous- 
tic analog of ZB in a macroscopic two-dimensional sonic 
crystal was observed. A general theory of ZB of a multi- 
band Hamiltonian is studied by David and Cserti 29 . Re- 
cently, Vaseghi^S et al. studied the effect of an external 
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perpendicular magnetic field on the ZB for both quantum 
wire and quantum dot within the two-sub-band model. 

The dimensionless parameter po — k$d dictates 
whether the motion will be oscillatory or not. It was 
shown in Ref. [14] that the motion will be oscillatory 
when po 3> I . It was also shown that the ZB vanishes 
when the strengths of the RSOI and DSOI are equal. 
This is due to an additional conserved quantity. 

In this work we study effect of the in-plane magnetic 
field on the ZB of electrons in a semiconductor QW and 
in a QD with the Rashba and Dresselhaus SOI. We ob- 
tain an analytical expression of the time-evolution of 
the position vector of an electron in a QW by using 
the Schrodinger picture. The advantage of using the 
Schordinger picture is to see the direct relation between 
the amplitude of the oscillatory motion and the Berry 
connection in momentum space. We find that in pres- 
ence of the magnetic field the ZB in a quantum well does 
not vanish even when the strengths of the Rashba and 
Dresselhaus spin-orbit interactions are equal. The ZB in 
a QW sustains even at very low value of kod due to the 
presence of the magnetic field. The trembling motion of 
an electron in a semiconductor QW with high Lande g- 
factor (e.g. InSb) sustains over a long time, even at very 
low value of kod. The time-period of the oscillation de- 
creases as the magnetic field strength is increased. Next, 
we study the ZB of an electron in a QD within the two 
sub-band model numerically The trembling motion does 
not fade away in time even when the magnetic field is ab- 
sent as well as when the strength of the SOIs is equal. 

This paper is organized as follows. In section II, we 
review the Hamiltonian and eigenstates of a 2DEG with 
both type of SOI in the presence of an in-plane mag- 
netic field. In section III, we derive time-evolution of the 
electron' position vector in the Schordinger picture. The 
numerical results and discussions are given in section V. 
Section VI contains the calculation and result of ZB of 
electrons in a GaAs/AlGaAs QD. The conclusion is pre- 
sented in Section VII. 



II. TWO-DIMENSIONAL ELECTRON GAS 
WITH SPIN-ORBIT INTERACTIONS 

We consider a 2DEG with the Rashba and Dresselhaus 
spin-orbit interactions in presence of an in-plane mag- 
netic field B = xB x + y By . The single-particle Hamilto- 
nian of this system is given by 

= 2^ li ( axPy ~ ° vVx > h ^ xVx _ ayPv > 

+ -^^ B { B X^ X +B V C7y), (1) 

where p and m* is the momentum and the effective 
mass of an electron, respectively. Here a and (3 are the 
strengths of the Rashba and Dresselhaus spin-orbit inter- 
action. The last term in the Hamiltonian is the Zeeman 
term due to the application of the in-plane magnetic field 



B, g* is the effective Lande-g factor and \x B is the Bohr 
magneton. The energy eigenvalues and the correspond- 
ing eigenstates of the syste m 31 ' 32 are, respectively, given 
by 

e± (k) = y ^ d ± Vc^TW*, (2) 

where C = ak y + /3k x + g* fi B B x /2, D = ak x + j3k v — 
g* [i B B y /2, and the spin eigenstates are 

!*- + > = 7lU-"')' '^> = 7i ("V) >< 3 > 

with 6 k = tan _x (a/D). 

In absence of the magnetic field, a new conserved quan- 
tity £ = <t x — <jy exists for a = / 3 33 i 34 . In this situation, 
the spin eigenstates \k,±) become independent of the 
wave vector k and therefore the ZB does not exist due 
to the absence of spin randomization 16,17 . The situation 
is completely different when the in-plane magnetic field 
is present. It is clear from Eq. © that the spin states 
are always function of k even at a = j3 as long as B ^ 0. 
Therefore, spin randomization occurs and we would ex- 
pect to see the ZB which is shown in the next section. 

A. Time-evolution of the wave packet and the 
Zitterbewegung 

We shall use the Schordinger picture to analyze the 
time-evolution of the electron position vector. We repre- 
sent the initial wave function of an electron by a Gaussian 
wave packet with initial spin polarization along the z axis 
as given by 

V(r,0)=-L J d 2 fca(k,0K k - r ( J) , (4) 

where a(k,0) = d/(V^)e~^ 2 ( k ~ k °) . Here, d and k 
is the initial width and the initial wave vector of the 
wave packet, respectively. The time-evolution of the ini- 
tial wave packet in the Schordinger representation can 
be obtained in the usual manner as tp(r, t) = U (t)tp(r, 0), 
where U (t) — e~ %m l n is the time-evolution operator. Af- 
ter doing some straightforward algebra we obtain, 

i)(r,t) = 1- [ d 2 fca(k,0)e lkr e- l 5^- t 

27T J 

x { cos[w(k)t] f J ) - eWk sia[w(k)t] ( J j }, 

(5) 

where w(k) = [E+Qs.) - J5_(k)]/(2/i). 

Equation ([5]) is nothing but the Fourier transformation 
of the following function: 

j.n j-\ ft -i-^rt ( COsL>(k)t] \ 

*(M)=a(k,0)e ^^-e^-sinKk^J- 
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The expectation value of position operator is then simply and 
given by 

(r(t)) =i /d 2 fc0t( k ,i)V k 0(k,i). (6) (y(t)) = 



Now it is straightforward to show that 



777, 



1 SLgLfrBs+By) 

2 (koy + ^) 2 + (Vkoy - 



x |l-cos[w(k)t]|. (7) 



The last oscillatory term cos[w(k)t] in Eq. ([7]) indicates 
the ZB. Two important observations can be made by ana- 
lyzing Eq. (J7J. First, the amplitude of the ZB is directly Ui Ms _ d(x 
proportional to the Berry connection (fc,±|z^j 
Vk^fe- Second, for a = (3 the amplitude (or the Berry 



We also calculate the expectation values of the velocity 
of electron in x and y direction. They are given by 



v x (t)} 



Vkffc- second, tor a = a tnc amplitude (or tne tserry „ , . r*n 7- 00 

• \r! .„ . , . , , _ V R -d 2 kl \ iii t -q 2 +2dqk 0y sin0 

connection) of the oscillation does not vanish when the — ^ e / a< P a 1 e 

in-plane magnetic field is present. But it vanishes if the " 
magnetic field is absent 



2/1 _ „2 



q 2 {l 



rf)sm<t> + q[e x +r)e y ) . (2aQ\ 

— sin 1 , 

Q V fid J' 



After substituting the expressions for a(k,0), V k 6» fc in Q 
Eq. ((7|) and considering the initial velocity of the wave (12) 
packet is along y-direction (i.e. k Qx — 0, k Qy — k ) one 
can easily obtain the following expressions and 



(x(t)) 



-d^kf. 



d</> dqe 
^ Jo Jo 

q 2 (l - rj 2 ) sin (j) + q(e x + rje y ) 



(2a 
1 — cos ( —tQ 



q -\-2dqkoy sin <fi 



KM 



(8) 



and 



{»(*)> 



fifci 



0y 



d 



-d-'kf. 



dqe~ q +2d i k °y sir "' i 
q 2 (r] 2 - 1) cos <f) + q(r)e x + e y ) 



27T JO -/(I 



(2a „ 
l-cos(-iQ 



(9) 



where Q 2 = (q cos <fi + ijq sin <fi — £ y ) 2 + (q sin tfi + ijq cos <fi+ 
e x ) 2 , T) = /3/a, q = kd, e x = g* (i B B x d/(2a) and e y = 
g*(L B B y d/(2a). 

From Eqs. ([5J and (J9j) one can infer that ZB is absent 
in the limit d — > 0. On the other hand, as d — > oo the 
Gaussian approaches to a delta function i.e. a(k, 0) = 
<5(k— ko) and in this limit we find an analytic expression 
for the ZB as given by 



(x(t)) 



1 

2(fc 0y 



(l-V 2 )koy + ^(B x 



VBy 



-) 2 + (V k 0y 



il - cos[w(fcoj,)i] j, 



(10) 



(Vy(t)) = 



hk, 



Of/ 



2vr 



2tt />oo 



— q -\-2dqkQy sin 



5 2 (t7 2 - 1) cos(j) + q{j]£ x + e y ) . (2aQ 

X — Sill 



Q 



V hd )■ 
(13) 



Here, v R = fid = 2a/ 'fi is the velocity corresponding to 
the RSOI. The corresponding currents are simply given 
by CM*)) = e(v x (t)) and (j y (t)) = e («„(*)), where e is 
the electronic charge. Eqs. ([9]) and (|13|) tell us that the 
ZB along y direction vanishes when in-plane magnetic 
field and the DSOI are absent simultaneously^ because 
in this case J Q V e 2dqk ° y COS( ^ sin (j)d<j) = 0. But when there 
is a finite in-plane magnetic field present the ZB in y 
direction does not vanish even at j3 = 0. So this is an 
effect of the in-plane magnetic field on ZB. 



B. Numerical Results and Discussion 

In this sub-section we evaluate time-evolution of 
the observables position vector and current density for 
different values of the parameters like magnetic field, 
kod, P etc. 

GaAs/AlGaAs QW: We consider GaAs/AlGaAs 
quantum well for which the effective Lande <?-factor is 
g* = —0.44. The value of the Rashba coefficient is taken 
to be a — 1.0 x 10 -11 eV-m. We set the condition dko = 5 
in all the cases. To investigate the time dependence of the 
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expectation values of position and current of electron, we 
numerically evaluate Eqs. ©, ©, {32) and CGI]). Here it 
should be mentioned that Eqs. ((9]) and (fTTj) contain two 
parts: the first depends linearly on time and the second 
one is oscillatory in time (responsible for ZB). The mag- 
nitude of the first term is very large compared to that of 
the second one. So if we plot this as a function of time we 
will get a straight line due to the dominating first term. 

We consider three cases corresponding to three differ- 
ent values of the parameter (3 namely (3 = 0, (3 = 0.5a 
and (3 = a. For each case we plot (x(t))/d, {jx(t))/eVR 
and (j y (t) — jo)/svR as a function of fit for different val- 
ues of magnetic field in Figs. [1-3]. Here we define the 
quantity jo a s jo = ehkoy/m*. From Fig. [1(a)], one 
can see that the amplitude of the ZB increases as the 
magnetic field increases from its zero value. It is also 
noticeable that the ZB pattern is more oscillatory with 
increasing magnetic field. The current shows (in Figs. 
[1(b), 1(c)]) the similar behavior as the position but it os- 
cillates about zero. From Figs. [2] and [3] we can see 
that as we increase the value of /3, the amplitude of ZB 
decreases. One important point is to be noted here that 
there is a definite phase difference between the currents 
in x and y direction when j3 =^ a as evident from Figs. [1] 
and [2]. But the situation is different when /3 = a. It can 
be seen from Fig. [3(b), 3(c)] the currents are oscillating 
in the same phase and this can be easily understood by 
analyzing Eqs. [121 113] . 
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FIG. 1: Here {x(t))/d, {jx{t))/ev R and (jy(i) — jo)/ev R are 
plotted as a function of Qt. In the all the cases we set /3 = 
for GaAs/AlGaAs QW. Here, solid line: (B x = 0, B v = 0), 
dotted line: (B x = I/\/2 T, B y = 1/V2 T), dashed line: 
(B x = 6 T, B y = 8 T). 



In Rcf. [14] it was shown that the ZB vanishes if a — (3 
in the absence of any external magnetic field. But Fig. 
[3] shows that the ZB is still present when a — (3. This is 
due to the application of the external in-plane magnetic 
field. 

InSb QW: Here, we consider InSb QW for which the 
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FIG. 2: Here (x(t))/d, {j x (t))/ev R and (j y (t) - jo)/ev R are 
plotted as a function of Qt. In all the cases we set /3 = 0.5a 
for GaAs/AlGaAs QW. Here, solid line: (B x = 0, B y = 0), 
dotted line: (B x = l/y/2 T, B y = l/y/2 T), dashed line: 
(B x = 6 T, B y = 8T). 
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FIG. 3: Here (x(t))/d, (j x {t))/ev R and (j y (t) - jo)/ev R are 
plotted as a function of Qt. In all the cases we set j3 — 0.5a 
for GaAs/AlGaAs QW. Here, solid line: (B x = 1/V5 T, B y = 
1/V2 T), dotted line: {B x = 3 T, B y = 4 T), dashed line: 
(B x = 6 T, By = 8 T). 



the effective Lande g-factor is very high (e.g g* = —50) 
as compared to GaAs/AlGaAs. The RSOI strength is 
taken to be a = 0.9 x 10~ n eV-m. We set here j3 = 0.5a 
and dko — 5. In Fig. [4] (x(t))/d is plotted with respect 
to ilt for different values of the magnetic field. But the 
situation is different from the GaAs/AlGaAs QW case. 
Although the amplitude decreases but the number of os- 
cillations contained in ZB within the same time range is 
quite large as we increase the magnetic field. Since the 
magnitude of g* is large the coupling between electron's 
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FIG. 4: (x(t))/d is plotted as a function of Sit for InSb QW. 
Here, we set /3 = 0.5a. 



This Hamiltonian can be re-written as 
H = hjj [a x a x + Oydy + ij 



—U 

V21V 

V21V 
9* 



'V 



a l) a y 



(B x a :r 



ByOy). 



(15) 



We consider only two lowest occupied energy states 
(ground state and first excited state) of a two- 
dimensional harmonic oscillator potential. This approx- 
imation is known as the "two sub-band model" . Within 
this approximation the Hilbert space spanned by the fol- 
lowing six basis vectors: |o,o,t), 10,0,4.), |i,o,t), |i,o,4-), 

|0, 1, t) , |0, 1, 4-)- Here, f and 4- represent the ^-component 
of the electron's spin vector. 

Within six basis vectors one can write the Hamiltonian 
in a matrix form as 



spin and magnetic field is strong. 

One can obtain more oscillation in ZB by increasing the 
magnitude of the parameter dko and we have shown that 
when dko » 1 the pattern is completely oscillatory as 
evident from Eq.(10). In all these cases it is observed that 
the ZB is transient in nature i.e it's amplitude decreases 
with time and it is a direct consequence of Lock's^ pre- 
diction that the ZB of electron will not be persistent in 
time if it is represented by a wave packet. 



III. ZITTERBEWEGUNG IN A QUANTUM 
DOT 



In this section we would like to study ZB of elec- 
trons in a semiconductor QD. We consider a 2DEG 
confined by an isotropic harmonic oscillator potential 
V(x,y) = (1 / 2)moj 2 (x 2 +y 2 ). In this context our Hamil- 
tonian reads as 



VyPx) 



t{<J x Px - cryPy) + —^ B (B x a x + By<7y). (14) 



We introduce the conventional harmonic oscillator 
creation and annhilation operators as a x — (x/l + 
ilPx/K)/V2, <4 = ( x /l - ilp x /h)/y/2, a y = (y/l + 
ilp y /H)/V2, and = (y/l — ilp y /H)/V2, where I = 
\J h/ (uiu) is the harmonic oscillator length. 
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The matrix elements are as follows: Hu = H22 
so, #33 = H44 = H 55 — H 66 — 2ftw — £1, H12 ■ 
#fo = H. 



#34 
H 



41 - 
#61 



#56 

-#i 

— 1 1 ' 

— ll^c 



'43 — "65 

— H32 
1 ; ■ _ 



: fkjj 

g*fi B (B x ~iB y )/2, H u 
(a-iP)/(lV2), #25 = #5*2 
(p + ia)/(lV2). 



Here £n = fow is the zero-point energy and £1 = 2huj 
is the 1st excited state energy of the two-dimensional 
harmonic oscillator. 

We want to determine the expectation value of the 
time-dependent position operator in this system. Let us 
consider at t = the system is in the ground state and the 
spin is oriented along the positive z-direction. This initial 
state is given by IV'(O)) = |0, 0,f), so the expectation 
value of the position operator is given by 

(x H (t)) = (Tp(o)\x H (t)\m) 

= (^(0)\e lHt/h x(0)e-' Ht/h \^(0)) 

= (0\VV^e lHt/h VV^xVV^e- lHt/h VV^\0) 

= (0\VV^e im ^VXVh- iHt / H VV^\0) 

= E(X>i^ ll ^)v4, (16) 
1=1 j=l 

where V is the diagonalization matrix which diagonalizes 
the Hamiltonian H and VV^ = I. Also, cjy = (Xi — Xj)/H 
is the beating frequency. 

The V matrix also diagonalizes U(t) and becomes 
U diag (t) = VU(t)V^ = e- iXit / h Iij and X = VxVl The 
position operator is given by x = A (a + a) )/y/2. Within 



6 



the above mentioned basis this can be written in a matrix 
form as 



V2 



/0 1 0\ 

1 





1 
\0 1 0/ 




with frequencies corresponding to the all possible energy 
eigenvalue differences of the Hamiltonian. In this case, 
there are six non-degenerate eigenvalues and we have six 
values of the energy differences. The number of beating 
frequencies is six. 




FIG. 6: Plots of (x(t))/\ vs cut for GaAs/AlGaAs QD for 
different values of /3. The value of magnetic field is kept fixed 
to (B x ,B y ) — (0.1,0.1)T. Figures (a),(b) and (c) correspond 
to /3 = 0, P = 0.5a and /3 = a respectively. 



FIG. 5: Plots of (x{t))/X vs ut for GaAs/AlGaAs QD for 
different values of magnetic fields. Here, Fig. (a): (B x = 0.1 
T, B y = 0.1 T), Fig. (b): (B x = l/yf(2) T, B y = 1/^(2) 
T), Fig. (c): (B x = 6 T, B y = 8 T). We consider only RSOI 
i.e. 13 = 0. 



A. Numerical Results and Discussion 

The ZB of a GaAs/AlGaAs quantum dot in the pres- 
ence of the in-plane magnetic field with SOI is investi- 
gated here. The value of the Rashba strength is taken as 
a = 1.0 x 10~ n cV-m and the zero-point energy of the 
harmonic oscillator potential is fixed to the value £o = 5 
meV. We find the expectation values of the position co- 
ordinate as a function of time which are plotted in Figs. 
[5] and [6]. In Fig. [5] we fix j3 = and vary the mag- 
netic field strengths. The magnetic field is kept constant 
and (3 is varied in Fig. [6]. The ZB in quantum dots is 
similar to the beating effect in the classical wave mechan- 
ics with different frequencies. The oscillatory motion is 
due to the superposition of individual oscillatory motions 



IV. CONCLUSION 

In this work we have investigated the effect of an ex- 
ternal in-plane magnetic field on the ZB of an electron 
in semiconductor QW and QD with Rashba and Dressel- 
haus spin-orbit interactions. For QW a general expres- 
sion of the expectation values of position coordinate and 
current due to ZB within the Gaussian wave packet is 
obtained. For QW case, the oscillatory quantum motion 
of electron which is represented by a wave packet shows 
transient behavior and this signature is a proof of Lock's 
argument. Another important point is that ZB does not 
vanish even at a = when a finite in-plane magnetic field 
is present. The y-component of current also performs ZB 
motion with finite magnetic field. We study the same 
problem for high Lande g-factor QW like InSb in compar- 
ison with low Lande g-factor QW like GaAs-AlGaAs. We 
have also studied the problem of ZB in a GaAs/AlGaAs 
QD numerically. The ZB in GaAs/AlGaAs QD is persis- 
tent in time. The ZB in quantum dots shows beating-like 
pattern and it is similar to the the beating effect in the 
classical wave mechanics. 
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